function [alpha,h,rmse_min]=ah_find(u0,T0,time,Temperature)
    s=0;i=1;N=100000;
    while(s<N)
        x=rand(1,2);
        a0=x(1);b0=x(2);
        rmse(i)=FittingError(a0,b0,u0,T0,time,Temperature);
        mine(i)=min(rmse);
        if(i==1)
            alpha=a0;h=b0;rmse_min=mine(1);
        else
            if(mine(i)<mine(i-1))
                alpha=a0;h=b0;s=0;rmse_min=mine(i);
            else
                s=s+1;
            end
        end
        i=i+1;
    end
end

